Accepted for publications in The Astrophysical Journal 

Preprint typeset using 1^1^^ style cmulatcapj v. 08/13/06 



0^ 
O 
O 



in 

< 
O 

6 



> 

00 

o 



m 
o 

o 



X 



A CHANDRA X-RAY ANALYSIS OF ABELL 1664: COOLING, FEEDBACK AND STAR FORMATION IN THE 

CENTRAL CLUSTER GALAXY 

C. C. Kirkpatrick\ B. R. McNAMARA^'^'•^ D. A. Rafferty"*, p. E. J. Nulsen'\ L. Birzan^, F. Kazemzadeh^ M. W. 

Wise"', M. Gitti'' and K. W. Cavagnolo^ 

Accepted for publications in The Astrophysical Journal 

ABSTRACT 

The brightest cluster galaxy (BCG) in the Abell 1664 cluster is unusually blue and is forming stars 
at a rate of ~ 23 Mq yr~^. The BCG is located within 5 kpc of the X-ray peak, where the cooling 
time of 3.5x10® yr and entropy of 10.4 keV cm^ are consistent with other star-forming BCGs in 
cooling flow clusters. The center of A1664 has an elongated, "bar-like" X-ray structure whose mass 
is comparable to the mass of molecular hydrogen, ^ 10^° Af© in the BCG. We show that this gas is 
unlikely to have been stripped from interloping galaxies. The cooling rate in this region is roughly 
consistent with the star formation rate, suggesting that the hot gas is condensing onto the BCG. We 
use the scaling relations of Birzan et al. (2008) to show that the AGN is underpowered compared to the 
central X-ray cooling luminosity by roughly a factor of three. We suggest that A1664 is experiencing 
rapid cooling and star formation during a low-state of an AGN feedback cycle that regulates the rates 
of cooling and star formation. Modeling the emission as a single temperature plasma, we find that 
the metallicity peaks 100 kpc from the X-ray center, resulting in a central metallicity dip. However, 
a multi-temperature cooling flow model improves the fit to the X-ray emission and is able to recover 
the expected, centrally-peaked metallicity profile. 

Subject headings: cooling flows — galaxies: clusters: individual: A1664 — galaxies: starburst — 
X-rays: galaxies: clusters 



1. INTRODUCTION 

The fate of hundreds to thousands of solar masses per 
year of gas thought to be condensing out of the X-ray at- 
mospheres of galaxy clusters is a problem that has puz- 
zled astronomers for more than three decades (Fabian 
1994). New Chandra and XMM-Newton observations 
have recast the problem in terms that have significant 
consequences for understanding the formation of galax- 
ies and super-massive black holes. First, high resolu- 
tion spectra of the thermal emission from cooling flows 
failed to show the characteristic recombination lines from 
highly ionized metals in the condensing gas (Peterson et 
al. 2003; see Peterson & Fabian 2006, for a review). The 
XMM-Newton observations did not rule out cooling al- 
together (eg., Sanders et al. 2008), but they were the 
first to show convincingly that the condensation rates 
must lie far below the levels predicted by pure cooling 
models. Second, images of cluster atmospheres revealed 
large-scale cavities and shock fronts associated with pow- 
erful AGN outbursts in central dominant galaxies (see 
McNamara & Nulsen 2007, for a review). Measurements 
of cavity sizes and their surrounding pressures provide a 
convenient and reliable gauge of the mechanical pV work 
(energy) expended by radio jets as they inflated cavities 
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against the surrounding gas pressure (McNamara et al. 
2000). Assuming the cavities are driven to their current 
locations primarily by buoyant forces (eg., Churazov et 
al. 2002), the mean jet power is comparable on average 
to the power required to quench cooling in the cores of 
galaxies and clusters (Birzan et al. 2004; Dunn & Fabian 
2006; Rafferty et al. 2006). Apparently, the supermas- 
sive black holes located in the nuclei of BCGs combined 
with abundant fuel accreting onto them provide a nat- 
ural and powerful feedback mechanism that is energeti- 
cally able to maintain most (but not all) of the cooling 
gas above a few keV (Briiggen & Kaiser 2002; Reynolds 
et al. 2002; Dalla Vecchia et al. 2004; Ruszkowski et al. 
2004a,b; Briiggen et al. 2005; Heinz et al. 2006; Vernaleo 
& Reynolds 2006). Thus the failure of XMM-Newton to 
find strong cooling lines can be attributed largely to AGN 
feedback. How this jet power heats the gas and the ex- 
tent to which other heating mechanisms such as thermal 
conduction, galaxy mergers, etc. are aiding it is poorly 
understood, and remains an outstanding issue (Zakam- 
ska & Narayan 2003; Dolag et al. 2004; Voigt & Fabian 
2004; Dennis & Chandran 2005; Poole et al. 2008). 

Evidence for residual cooling can be inferred by the un- 
usually high star formation rates and reservoirs of cold 
gas found in BCGs. Although this star formation is fre- 
quently attributed to stripping from the occasional gas- 
rich galaxy (eg., Holtzman et al. 1996), several surveys 
undertaken over the last few decades have tied the pres- 
ence of nebular line emission, molecular gas, and star 
formation to cooling flows (Johnstone et al. 1987; Mc- 
Namara & O'Connell 1989; Crawford et al. 1995, 1999; 
Cardiel et al. 1998; Edge et al. 2002; Salome & Combes 
2003; Rafferty et al. 2006; Donahue et al. 2007; Edwards 
et al. 2007; Bildfell et al. 2008; O'Dea et al. 2008; Raf- 
ferty et al. 2008). Despite these ties, the enormous gap 
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between the cooling and star formation rates raised very 
serious issues. This situation again changed dramati- 
cally when XMM'Newton^s downward-revised condensa- 
tion rates were shown to be in near agreement with the 
star formation rates in many, but not all, systems (Mc- 
Namara 2002; Hicks & Mushotzky 2005; RafFerty et al. 
2006; O'Dea et al. 2008). 

Recently, evidence tying the cooling of hot halos to 
star formation was found by RafFerty et al. (2008), who 
showed that star formation ensues when the central cool- 
ing time falls below a remarkably sharp threshold value 
of ~ 5x 10* yr or equivalently a central entropy of 30 keV 
cm^. This and two other criteria: the X-ray and galaxy 
centroids lie within ~ 20 kpc of each other, and the X- 
ray cooling luminosity exceeds the jet (cavity) power, 
appear to govern the onset of star formation in cooling 
flows. Cavagnolo et al. (2008) found a similar threshold 
for the onset of Ha and radio emission (see also Hu et 
al. 1985). The Ha threshold is certainly related to but 
not identical to the star formation threshold, as we show 
in this paper. These new results are difficult to under- 
stand in the context of mergers or stripping, but instead 
they tie the presence of star formation and AGN activ- 
ity closely to cooling instabilities in the hot atmospheres 
(Rafferty et al. 2008; Voit et al. 2008; Soker 2008). 

These developments have broad implications for under- 
standing the formation and evolution of galaxies. Out- 
flows driven by gas accretion onto massive black holes 
may have regulated the growth of bulges giving rise to 
the observed correlation between bulge mass and super- 
massive black hole mass in nearby galaxies (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000; Haring & Rix 2004). 
In simulations, AGN feedback at late times in the so- 
called "radio mode," of which feedback in cooling flows 
is thought to be the archetype, is able to suppress cool- 
ing and star formation in the hot halos of giant ellip- 
tical galaxies and BCGs (Bower ct al. 2006; Croton et 
al. 2006). This process may explain the turnover at the 
bright end of the galaxy luminosity function and the lack 
of bright blue galaxies expected in standard ACDM mod- 
els. Cooling flows permit study of this process in detail. 

Abell 1664 is a good candidate for studying star for- 
mation and cooling in galaxy cluster cores. It is a cool- 
ing flow cluster (Allen et al. 1995) located at rcdsliift 
z = 0.128 (Allen et al. 1992) and has an Abell richness 
class of 2 (Abell et al. 1989). Based on ROSAT data, the 
X-ray luminosity in the 0.1-2.4 keV band is 4.1x10''^ erg 
s-i (Allen et ai. 1992). A1664 hosts a BCG known to 
be a very bright Ha emitter with a Imninosity of at least 
1.6x10^^ erg s~^, and forming stars at a rate of 14-23 
Mq yr-i (Wilman et al. 2006; O'Dea et al. 2008). In 
the X-ray waveband, the cluster is regular in shape on 
large scales. The core shows signs of a surface brightness 
break along the eastern edge at a radius of 30" from the 
center and an elongated surface brightness peak at the 
center. 

In this paper, we present our analysis of the recent 
Chandra observations and deep U- and i?-band images 
of the cluster and its core. Section 2 describes the data 
set and reduction process. In Section 3, we describe our 
analysis of the X-ray and optical data. In Section 4 we 
interpret the data. Finally, we summarize in Section 5. 
Throughout the paper we assume a ACDM cosmology 
with Ho = 70 km s'^ Mpc'S VLm = 0.3, and = 0.7. 



The angular scale is 2.29 kpc arcsec" at the redshift of 
the cluster. All uncertainties quoted are 90% confldence 
intervals. 

2. OBSERVATIONS & DATA REDUCTION 

Abell 1664 was observed for 37 ks on 2006 December 
4 (ObsID 7901) using the ACIS-S3 CCD on-board the 
Chandra X-ray Observatory. The focal plane tempera- 
ture during the observation was kept at -120°C. CIAO 
version 3.4 was used for all data reduction and prepara- 
tion with version 3.4 of the calibration database. Stan- 
dard filtering was applied to eliminate background flares, 
producing a loss of 0.4 ks from the total exposure time. 
Point sources have been identified and removed using 
wavedetect and confirmed by eye. Time-dependent gain 
and charge transfer inefficiency corrections have been ap- 
plied. Blank-sky background files used for background 
subtraction were normalized to the source image count 
rate in the 9.5 — 12 keV band. When extracting spectra 
from a region, a background spectrum is extracted from 
the equivalent region of the blank-sky background event 
file. 

Optical images of the cluster's BCG were obtained 
with the Kitt Peak National Observatory's 4 m telescope 
on 1995 February 2 using the T2KB CCD camera. The 
seeing at the time was approximately 1.3". Images were 
taken through U and Gunn R filters with total exposure 
times of 2400 s and 800 s, respectively. These filters were 
chosen to avoid contamination by emission lines from 
[OH] 3272 and Ha 6563 at the redshift of the cluster. 
Additional imaging was done using an i?-band filter that 
includes Ha emission with a total exposure time of 800 
s. The individual images have the bias level subtracted 
and have been flat-fielded using twilight sky exposures. 
Finally, the images for each band were combined into 
single U and R images and color corrected for use in our 
analysis. An Ha-|-[NII] emission map was created by tak- 
ing the scaled difference between the normalized R- and 
Gunn i?-band images. 

3. DATA ANALYSIS 
3.1. Images 

We present in Figure 1 the i?-band DSS image of 
A1664. Overlaid are the X-ray contours spanning the en- 
tire ACIS-S3 CCD (8' X 8'). On large scales, the X-ray 
emission is fairly regular, but with asymmetries to the 
south and east. The inner region is more complex. The 
cluster core image in the 0.4 - 7.0 keV band is presented 
in the upper left panel of Figure 2. A bar- like structure 
is shown as the dark region approximately 22.5" x 10" 
(52 X 23 kpc) in size associated with the BCG, marked 
A. To the east, a brightness edge, marked B, is seen at 
a radius of 30" from the cluster center. In the upper 
right panel of Figure 2, Ha emission is aligned with the 
"bar" region, though the bar is more extended than both 
the Ha and star forming region. The feature marked C 
is thought to be a disturbance created by a neighbor- 
ing galaxy (Wilman ct al. 2006). In the bottom left and 
right panels of Figure 2 we present the R- and [/-band 
images, respectively. The i?-band shows that the BCG's 
core docs not share the same position angle on the sky as 
the outer envelope. In the [/-band, there is more struc- 
ture. There are two distinct clumps north and south 
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Fig. 1. — DSS R-band image of A1664. The X-ray contours show 
the emission spanning the entire ACIS-S3 chip. North is toward 
the top; East is to the left. 

of the center. We will discuss possible origins of these 
features in the following sections. 



3.2. X-ray Properties 

In the following sub-sections, we present our analysis 
of the Chandra observation of A1664. Our purpose is to 
obtain the mean properties of the cluster including the 
gas temperature, density, metallicity, cooling rate, and 
how these properties vary with radius. 

3.2.1. Surface Brightness Analysis 

We determined the center, ellipticity and position an- 
gle of A1664 on the sky by fitting a two-dimensional 
Lorentzian surface (Wise et al. 2004). The centroid is 

found at R.A. = 13''03"42^5, DEC = -24°14'43".96. 
This centroid is offset from the cluster emission peak by 
3.6" to the NE. The average ellipticity was found to be 
0.22 with a position angle of 25°. Both are consistent 
with values found by Allen et al. (1995). When dividing 
the data into annular regions for measuring radial trends, 
the centroid, ellipticity, and postion angle derived here 
were used. 

An exposure map created using mkexpmap was divided 
into the image. The surface brightness was calculated by 
measuring the flux through annuli 5" in width. We fit 
the surface brightness profile with a beta model using 
the Sherpa model betald. A single beta model does not 
adequately fit the cluster surface brightness. The model 
adequately fits the outer regions of the profile, but not 
the elevated brightness of the core. Adding a second 
component to fit the cooling cusp improves the fit. The 
form for the double-beta model is 

/(r) =Ji(1 + 4)""'^+5+/2(1 + 4)"'''=+5+/b, (1) 

where and are the amplitudes and core radii. The 
constant compontent, Ib, was added to represent the 
background. All parameters from the fit can be found 



in Table 1. The surface brightness profile presented in 
Figure 3 includes the fit from the double-beta model. 
The smooth model does not fit the central region well, 
which is most likely due to the complicated structure of 
the core. 

The residual map formed by subtracting the double- 
beta model is presented in Figure 4. A spiral pattern is 
apparent near the center. Similar spiral patterns have 
been seen in other clusters (Sanders & Fabian 2002; 
Clarke et al. 2004). The spiral is commonly attributed to 
the angular momentum of stripped gas from galaxies or 
groups. Other explanations involve merging dark matter 
halos with trapped gas (Clarke et al. 2004; Dupke et al. 
2007). 

The spiral feature is common in systems with cold 
fronts (see Markevitch & Vikhlinin 2007, for a review). 

Most cold fronts are likely due to oscillations of the dark 
matter or gas in the cluster core (Tittley & Henriksen 
2005). This "sloshing" can be triggered by a tidal or 
pressure disturbance caused by an infalling subcluster 
or possibly a powerful AGN outburst. The sloshing gas 
can acquire angular momentum, which is why the spiral 
feature is associated with many cold fronts (Ascasibar 
& Markevitch 2006). We present evidence that A1664 
contains a cold front in section 3.2.5. 

3.2.2. Integrated Spectral Analysis 

For our spectral analysis, we begin by extracting the 
total spectrum of the cluster core using dmextract. The 
elliptical region used in the extraction was chosen with 
semi-major axis of 161" to stay within the edge of the 
ACIS-S3 CCD's footprint. Weighted response files were 
created for the extracted regions using mkacisrmf and 
mkwarf. Four different models were used to fit the prop- 
erties of the X-ray gas over the energy range 0.4 keV 
to 7 keV with the data grouped with a minimum of 20 
counts per bin. Each model type was fit twice, once 
with the column density (Nh) as a free parameter and 
once with column density set at the Galactic value of 
8.95x10^° cm-2 (Dickey & Lockman 1990). All mod- 
els allowed temperature, abundance, and normalization 
vary. The simplest model used was a single temperature 
(if) plasma with absorption (WABSxMEKAL). To test 
for multi-phase gas, a two temperature (2T) model was 
used by adding a second MEKAL component to the IT 
model while tying the abundance parameters together. 

We have also tested whether a cooling flow 
fits the spectrum by adding a cooling flow 
component to the single temperature model 
(WABSx(MEKAL-l-MKCFLOW)). The tempera- 
ture of the thermal component was tied to the upper 
temperature of the cooling flow component. Two 
separate fits were found by testing different scenarios 
based on the MKCFLOW model temperatures. For 
the first scenario, we allowed both upper and lower 
temperatures to vary (VC), as well as abundance and 
normalization. In the second case, we set the lower 
temperature, kTiow, to 0.1 keV (e.g., Wise et al. 2004), 
allowing the gas to cool fully (FC), i.e., to below the 
temperature range detectable by Chandra. From this 
we found an upper limit on the mass condensation rate 
for the cluster within in the ACIS-S3 chip to be 137 Mq 

The results are in Table 2. The fits found for each 
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Fig. 2. — Upper left: Inner 74"x74" region in the 0.4 - 7.0 keV band adaptively smooth using CIAO tool aconvolve. Upper right: Ha 
map of the same region. Lower left: R-band image of the same region in false color to show structural detail. Lower right: C/-band image 
of the same region in false color. The labels are as follow: Point A is the X-ray bar, point B is the surface brightness edge, and point C is 
Hcf disruption being caused by a neighboring galaxy. Points 1 and 2 are clusters members that may be interacting with the BCG. 



model range in reduced values from 1.1 to 1.5. Multi- 
ple temperature models are favored over the single tem- 
perature model, which is expected in a cooling-flow clus- 
ter. The temperature variations will be examined more 
closely in section 3.2.5. All models with Nh allowed to 
vary proved better fits than the models set to the Galac- 
tic value, suggesting an error in the interpolation of the 
foreground value or error in the calculation of the buildup 
of hydrocarbons on the optical filter. The latter possi- 
bility is unlikely. A detailed CXC memo by Vikhlinin^ 
claims the calibration of the buildup is known to less 
than 5% at low energies. 

To understand how this absorption is distributed 
across the field of view, we created elliptical annuli con- 



taining about 5000 net counts. Fitting each of these re- 
gions with a single temperature model over the same en- 
ergy range as before, we produced the Nh radial profile 
shown in Figure 5. The dotted line indicates the Galactic 
value for hydrogen column density. The profile shows no 
radial trend, but overall the column density is above the 
Galactic value in all regions. The average column den- 
sity for this region is better fit with the value 1.22x10^"'^ 
cm~^. This value is in agreement with the fitted val- 
ues for the integrated spectrum. For all models used in 
this analysis, a significantly better fit is found using this 
average column density rather than the Galactic value. 
Therefore, we have adopted it for the remainder of this 
paper. 



See http://hea-www. harvard. edu/'^alexey/acis/mcmos/. 



3.2.3. The X-ray Bar 
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TABLE 1 
Beta-Model Fit 



ri 7-2 
(arcsec) f3i /i^- (arcsec) 02 Ib"^ xlorm d.o.f. 

89.83l||^ 0.78tg:g^ 75.7t\l 23.65lgJ^ O-Sltaoi 12201^? 6.93tg;g^ 2.57 228.48 89 



^Units of 10-9 phot sec ^ cm ^ arcsec ^ 

The bar region, indicated in Figure 2 (A), stands out 
as a strong peak in surface brightness. We determined 
its dimensions from its average surface brightness pro- 
files in the directions perpendicular and parallel to its 
major axis. Using the full-width at half-maximum of the 
average peaks found in Figure 6, the length is 22.5" and 
width is 10". These dimensions were used to isolate its 
spectrum. A total of 3350 net counts were detected in 
the extracted region. The center of a cluster suffers the 
most from projection effects, so contributions from four 
annular regions outside of the bar containing 15,000 net 
counts each were fitted for deprojection. The resulting fit 
of the PROJCT(WABSxMEKAL) model for the bar is 
presented in Table 3. The electron density of the region is 
0.066 cm^'^, corresponding to a total mass of 3.21x10"'^'^ 
Mq. If the gas in the bar was all cooling to low temper- 
atures on its cooling timescale, it would be depositing 
about 91 Mq yr"^. The density, pressure, entropy, and 
cooling time are all consistent with the values fit for the 
central region of the radial analysis, found later in section 
3.2.5. 

We will now consider the possibility that the bar is 
actually shaped more like an edge-on disk. Before we 
assumed that the bar is prolate, with its major axis in 
the plane of the sky. The edge-on disk is oblate, with 
its minor axis assumed to lie in the plane of the sky. 
With a larger volume because of the oblate geometry, 
a lower density will be measured. These differences are 
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Fig. 3. — Radial surface brightness profile in the 0.4 - 7.0 keV 
band. The fit shown represents the two components of the double- 
beta model and the best fit background value. The panel below 
shows the residual deviation from the fit. 



presented in Table 3. With these assumptions the density 
and pressure decrease to 0.0443 cm"^ and 2.42x10"^° 
erg cm"^, respectively. The entropy and cooling time 
rise to 13.6 keV cm^ and 5.3x10^ yr, respectively. The 
mass increases to 7.22x10^" Mq, more than double the 
mass contained within the bar geometry. The assumed 
geometry does not affect our conclusions. 

The CO-derived molecular gas mass for A1664 is 
4.4x10^° Mq (Edge 2001). The total mass of the X- 
ray emitting gas in the bar is comparable to the cold 
molecular gas that resides at the same location. A study 
by Wilman et al. (2006) has shown that the CO and Ha 
emission are tracing the same gas, which is shown in the 
upper right panel of Figure 2. The X-ray bar structure 
is associated with this same disturbed region. 

3.2.4. Metal Abundances 

The metallicity analysis was carried out by creating 
adaptivcly sized regions containing an average 5000 net 
counts out to the edge of the CCD (161"). AU re- 
gions have their own weighted response files created using 
mkacisrmf and mkwarf. After the spectrum for each re- 
gion was extracted, each was modeled as a single temper- 
ature plasma (MEKAL) with an absorption component 
(WABS). The model was fit over the energy range 0.4 to 
7.0 keV, with Nh set to our adopted value. Temperature, 
metallicity, and normalization were allowed to vary. The 
abundance ratios used in this analysis are from the most 
current photosphcric data of Grevesse & Sauval (1998). 




Fig. 4. — A 5'x5' residual map of the double-beta model sub- 
tracted image. The image was adaptively smooth using CIAO tool 
aconvolve before subtraction. 
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TABLE 2 
Integrated Spectral Fits 



Model 


kT 
(keV) 


(keV) 


z 


(1022 cm-2) 




x2/d.O./ 


IT 
IT 
2T 
2T 
VC 
VC 
FC 
FC 


,7+0.87 
3--'-' -0.32 
c 9^+0.34 

4 Q1+0.11 

o.oo_3 Q4 
4 23+° " 


1 69+"-" 

l.DO_g 22 

,7+0.04 

1 ,,+0.08 

J..OO_Q gg 
0.1 
0.1 


46+0 03 
0.54troi 

o.37i;!-;« 
o.49+i!i;i^ 
o.46t;!:ii 

n 47+0.03 

n r:fi+0.03 
<J.OD_() Q3 


U.iiD„Q gQ3 

0.0895 
120+"'004 

0.0895 
n loQ+o oos 

0.0895 
n 1,9+0.004 

U.10Z_Q gg4 

0.0895 


145111 
367121 

137+1^ 
29+10 


481.26/368 
564.01/369 
402.25/366 
503.93/367 
397.59/366 
503.19/367 
400.63/367 
556.21/368 



The projected metallicity profile is presented in the left 
panel of Figure 7 represented by the solid points. The 
abundance rises from the center of the cluster at 0.54 so- 
lar to 0.75 solar at 21" and remains constant out to about 
45". The metallicity then declines to approximately 1/3 
of the solar value at a radius of 100". The central dip in 
metallicity is unusual, although central abundance dips 
have been observed in several clusters such as Centau- 
rus A (Sanders & Fabian 2002), Perseus (Schmidt et al. 
2002), and M 87 (Bohringer et al. 2001). Usually the 
metallicity rises to the center in cooling flows (De Grandi 
& Molendi 2001). Our current understanding of the cen- 
tral metal enhancement of cooling flow clusters posits 
that SN la in the central galaxy are enriching the gas 
with iron and other metals. The peak metallicity is ex- 
pected to be centered on the BCG. These rare systems 
with central dips may be explained by resonance scatter- 
ing (Gilfanov et al. 1987), multiphase gas (Buote 2000), 
warm+hot ISM (Buote 2001), or bimodal metallicity dis- 
tributions (Morris & Fabian 2003). 

We explore the metallicity profile first by taking into 
account projection effects. A deprojected fit was carried 
out using 5 radial bins, containing roughly 10000 net 
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Fig. 5. — Hydrogen column density profile (projected). The 
dashed line represents the Galactic value. 



counts each. The spectrum for each region was modeled 
as a single temperature plasma with foreground absorp- 
tion, but now with the additional projected component. 
The open circles in Figure 7 representing the deprojected 
profile show that the abundance dip becomes more pro- 
nounced. Overall the 5 bins generally remain within er- 
rors of the projected profile. The peak at 45" lies at a 
value of 0.93 solar. To within errors, this value is consis- 
tent with the projected value at the same radius. 

We now explore the possibility of multiphase gas and 
low energy absorption due to warm gas. Molendi & 
Gastaldello (2001) have shown that a single temperature 
model can bias a measurement to lower mctallicities. It 
was previously demonstrated by Buote (2000) that this 
bias is due to the abundance measurement being sensi- 
tive to the Fe-L complex at low temperatures. In a multi- 
phase medium, the Fe-L complex is excited by a range of 
temperatures around 1 keV, resulting in a broader spec- 
tral shape. A single, average temperature would treat 
this broad peak as a lower metallicity feature, resulting 
in an underestimation of the actual abundance. 

Using the same regions and model we used for the pro- 
jected analysis, we added a cooling flow component to 
the model (WABSx(MEKAL-hMKCFLOW)) with the 
upper temperature tied to the temperature of the ther- 
mal component. The lower temperature was set at 0.1 
keV to allow full cooling. A significant improvement in 
the fit was found for the inner two regions when adding 
this component. An F-test shows that the outer six re- 
gions improve slightly or remain the same. The reduced 
chi-squared values for both models are presented in Table 
4. 

The new abundance profile with the inner two re- 
gions modeled with cooling flows is shown in the right 
panel of Figure 7. The central metallicity dip has largely 
gone away. To within the errors, the profile is centrally 
peaked. The upper limit on the mass condensation rate 
for these two regions was found to be 56±10 Mq yr~^. 
In an attempt to improve the fit of the inner most region 
further, we assume that there may be reservoir of warm 
gas (~ 10^ K), being deposited from the cooling flow, act- 
ing as a source of extra absorption in the lower energy 
range of the spectrum (0.4 - 0.8 keV). There is a known 
reservoir of cold gas that is comparable in mass to the hot 
X-ray emitting gas associated with the BCG (discussed 
in sec 3.3.1). Our assumption is that if the cooling flow is 
directly feeding the cold gas, then there would be gas in 
the intervening temperature range, though the amount 
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Fig. 6. — Left: Average surface brightness (in total counts) parallel to the length of the bar. Right: Average surface brightness (in total 
counts) parallel to the width of the bar. The red lines indicates the profile directly through the bar. The blue line is the average profile of 
the region immediately outside of the bar for comparison. 



TABLE 3 
Spectral Fits for Core Region 



Geometry 


(10-2 cm-3) 


P 

(10-10 erg cm-3) 


5 

(keV cm2) 


^cool 

(108 yrs) 


M 

(IQiO Mq) 


Mx 


Bar 
Disk 


6 64+°-^'' 
'-'•°^-0.28 

4 43+0-17 

^•^'-'-0.19 


O.DO_Q 23 

2 42+° " 


10 4+°-^ 
13 6+° '' 


n c + 0.4 
''•^-0.4 
c: Q+0.6 
°"^-0.6 


o 91+0.12 
■^■^^-0.14 
•J i-\(-j-(-0. 28 
'•^^-0.30 





expected is small due to extremely short cooling times 
and therefore the effect should be minimal. This extra 
absorption was represented using the EDGE component 
in XSPEC. A visual inspection of the spectrum at lower 
energies shows the fit is no longer systematically higher 
than the data, which increases the metallicity measure- 
ment to the solar value. However, by introducing two 
new free parameters to the model, the chi-squared value 
did not improve significantly. 

3.2.5. Physical State of the Gas 

A radial temperature analysis was performed by ex- 
tracting regions from the same area of the CCD as the 
metallicity analysis. The sizes of the annuli were adjusted 
so that each contained approximately 2000 net counts. 
The radial bins were fit in projection following the same 
procedure and using the same model as described above 
in the metallicity analysis. 

Shown in Figure 8 as asterisks, the temperature reaches 
a minimum of 2.2 keV at the center of the cluster. At a 
radius of about 18", the temperature rises from 2.6 ±0.1 
keV to 2.9 ± 0.2 keV and drops back down to 2.5 ± 0.2 
keV in the next annulus. This feature may be due to an 
asymmetric temperature structure caused by the eastern 
cold front. From there, the profile rises steeply to 4.6 
keV at about 45" and reaches its peak temperature of 
5.7 keV at around 100". We note that when fitting the 
temperatures with the abundances tied to those derived 



using the cooling flow model in section 3.2.4, the changes 
in the best fit values are imperceptible. 

The deprojected temperature profile was divided into 
8 elliptical annuli. The temperature was determined us- 
ing the PROJCT component in XSPEC, as before. The 
open circles represent the deprojected temperature pro- 
file in Figure 8. Apart from the central bins, the de- 
projected profile agrees with the projected profile. The 
central temperature decreases to 1.7 keV, making it 0.5 
keV cooler than found in the projected profile. Between 
the second and third radial bin, the temperature jump is 
still apparent in the temperature structure. 

Density, pressure, entropy, and cooling-time profiles 
have been derived from the same radial bins. The den- 
sity profile, obtained from the deprojection of the sur- 
face brightness profile, is shown in the upper left panel 
of Figure 9. The core reaches an electron density of 0.07 
cm--^, which is comparable to that found in the core of 
the Perseus cluster (Fabian et al. 2006). The profile de- 
creases radially to below 0.002 cm^'' in the last annulus 
at a radius of 160". 

The pressure profile was calculated as 2kTne- The 
upper right panel of Figure 9 shows the pressure peaking 
at the center at a value of 5xl0~^° erg cm"'^. At 20", 
there is a break in the profile. This is related to the 
temperature decrease described above. 

The entropy profile in the lower left panel of Fig- 
ure 9 was calculated using the standard relationship 
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Fig. 7. — Left: Projected (asterisks) and deprojected (circles) metallicity profile of the ICM in solar units using a single temperature 
model. An unusual dip occurs in the central ~ 50" . Right: Metallicity profile of ICM in solar units. The inner two points were derived by 
including a cooling flow model to obtain a better metallicity estimation. 



TABLE 4 

Reduced Chi-squared Fits for Radial Metallicity Profile 









Reduced chi-squared for region (d.o.f.) 






Model 


Reg. 1 


Reg. 2 


Reg. 3 Reg. 4 Reg. 5 Reg. 6 


Reg. 7 


Reg. 8 


IT 
CF 


1.15(133) 
0.98(132) 


1.35(146) 
1.07(145) 


1.06(148) 1.20(153) 0.99(156) 1.17(152) 
0.97(147) 1.15(152) 0.92(155) 1.17(151) 


0.95(206) 
0.93(205) 


1.12(234) 
1.10(233) 



S = kTn,. . The profile shows evidence for feedback 
in the core. A cluster without heating would have an en- 
tropy profile that decreases roughly as a power law to the 
center (Voit et al. 2002; Kay et al. 2004). We find that 
the profile decreases inward and then flattens at small 
radii around 13 keV cni^. This entropy floor is normally 
found in cooling flow clusters (Voit & Donahue 2005) 
where AGN activity is thought to be flattening out the 
profile by periodically releasing energy back into the sur- 
rounding gas. We present the cooling time as a function 
of radius in the lower right panel of Figure 9. The short 
central cooling times are characteristic of a cooling flow 
cluster. The decreasing cooling time towards the center 
reaches its minimum value of 3.6x10* yr, which is con- 
sistent with star forming cooling flows. The sample of 
Raffcrty et al. (2008) showed that star formation is pref- 
erentially found in systems with a central cooling time 
less than approximately 5x 10* yr, and with entropy less 
than 30 keV cm^. 

We pointed out in section 3.2.1 that the central struc- 
ture is asymmetric. A cold front could explain the struc- 
ture in the global temperature profile. We have analyzed 
the eastern half of the cluster, between 99 and 273 de- 
grees from west on the sky, towards the brightness edge 
pointed out in Figure 2 (B). The surface brightness pro- 
file in the direction of this region is shown in the up- 
per left panel of Figure 10. The break in the surface 



brightness profile at approximately 30" corresponds to a 
density and temperature discontinuity shown in the up- 
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Fig. 8. — Projected (asterisks) and deprojected (circles) X-ray 
temperature of the ICM. 
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Fig. 9. — Upper left: Density profile of tfie ICM. Upper right: Pressure profile of the ICM. Lower left: Entropy profile of the ICM. Lower 
right: Cooling time profile of the ICM. Bins below the dotted line (7.7x10^ yr) are considered to be within the cooling region. 



per right and lower left panel of the Figure 10. Aeross 
the front, away from the center, we find that density falls 
from 0.018 cm~^ to 0.008 cm~^. The temperature jumps 
at the front from about 3 keV to about 5.5 keV. How- 
ever, the pressure profile in the last panel of Figure 10 is 
continuous and shows no break. This is characteristic of 
a cold front, possibly induced by sloshing of the ICM in 
the core due to a recent merger. A1664 hosts a radio relic 
in the outer halo of the cluster (Giovannini et al. 1999; 
Govoni et al. 2001). Radio relics are thought to be a sig- 
nature of mergers (Enfilin et al. 1998; Enfilin & Briiggen 
2002), strengthening the case for merger induced sloshing 
in A1664. 

3.3. Star Form,ation in the BCG 

The visual light of A1664's BCG, shown in Fig- 
ure 2, is centered at R.A. = 13''03'"42«.54, DEC = 



-24 14'42".05, which hes 2.18" NE from the X-ray cen- 
troid. The R- and C/-band surface brightness profiles 
are presented in the left panel of Figure 11. We con- 
structed these profiles using the ellipse routine in IRAF 
with elliptical annuli of ellipticity 0.46 and position an- 

o 

gle 10 from north on the sky, based on the best fit 
to the BCG envelope. Nearby galaxies projected onto 
the envelope were masked out to avoid light contami- 
nation. The images were flux calibrated using Landolt 
standards. The Galactic foreground extinction correc- 
tion was calculated following Cardelli et al. (1989), as- 
suming a foreground color excess of E(i? — t^) = 0.141. 
Evolution and K-corrections from Poggianti (1997) were 
applied in both filters as -0.143 and 0.130 magnitude in R 
and -0.161 and 0.328 magnitude in [/, respectively. The 
error bars represent the statistical errors and the dashed 
lines represent the systematic error confidence intervals. 
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Fig. 10. — Upper left: X-ray surface brightness profile of the eastern half of A1664 in the 0.4 - 7.0 keV band. Upper right: Eastern half 
density profile of the ICM. The vertical line at ~ 30" represents where the surface brightness break is located. The density is discontinuous 
at this point as well. Lower left: Eastern half temperature profile of the ICM. The temperature sees a ~ 2 keV jump at the breaJf. Lower 
right: Eastern half pressure profile of the ICM. The pressure remains continuous across the surface brightness break. 



The method of determining these is described in McNa- 
mara & O'Connell (1992). 

In the right panel of Figure 11, we present the 
U—R color profile. The inner 10" shows a positive gra- 
dient (bluer toward the center) with a central rest-frame 
color of 1.1. The envelope of the BCG reddens to a color 
of approximately 2.3. The central color is about one mag- 
nitude bluer than a normal giant elliptical, which have 
colors of 2.3 - 2.6 (Peletier et al. 1990). This indicates 
star formation. 

The t/-band image, after subtraction by a smooth stel- 
lar model, is presented in the left panel of Figure 12. Two 
distinct peaks can be seen in this image north and south 
of the center. To the right of this in Figure 12, we show 
the color map on the same scale. The normal colors are 
in white and the unusually blue regions are blue. The 



blue features are concentrated to the center with trails 
of emission leading north, south, and east. The trails 
to the north and south appear to end at or near neigh- 
boring galaxies. The two galaxies labeled 1 and 2 in 
Figure 2 are both red and have an average color of 2.5 
and 2.2 respectively. Galaxy 1 has an absolute visual 
magnitude of approximately -22.7 {2Li,) and galaxy 2 is 
approximately -21.1 (O.SL^,). These appear to be giant 
ellipticals which are normally devoid of cold gas, so their 
relationship to the star formation in the BCG, if any, is 
ambiguous. They are unlikely donors of the 10^° Mq 
of cold gas feeding star formation in the BCG. 

4. DISCUSSION 
4.1. AGN Feedback 
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Fig. 11. — Left: R- (upper) and (7-band (lower) surface brightness profiles of the BCG. Statistical errors are imperceptible at most 
points. The dashed lines represent the systematic uncertainty associated with the sky background subtraction. Right: U — R color profile 
of the BCG illustrating the central ~ 10" starburst region. Statistical errors and systematic error envelope are shown. 




AGN periodically inject energy to heat the ICM and 
quench cooling of the surrounding medium (Birzan et al. 
2004; Dunn & Fabian 2006; Rafferty et al. 2006). To 
determine whether heating is dominating in A1664, we 
compare the mechanical power of the AGN to the cool- 
ing power of the ICM. AGN power can be evaluated as 
the pV work required to inflate a cavity against the sur- 
rounding pressure of the hot ICM. In A1664, there are 
no detected cavities in the X-ray. We can, however, esti- 
mate the AGN power using the scaling relations found by 
Birzan et al. (2008) between jet mechanical power and 
radio synchrotron power. These relations were derived 



using a sample of radio galaxies with X-ray cavities. All 
results referred to in this discussion are presented in Ta- 
ble 5. 

The radio flux emerging from the nucleus of the BCG 
was obtained from the NVSS catalogue. The source is 
unresolved in the 45" FWHM beam. Its 1.4 GHz flux 
36.4 mJy corresponds to a power of 1.57x10^^ W Hz^^. 
The bulk of the energy from FR I radio sources emerges 
in the form of mechanical energy. Using the 1.4 GHz 
relationship of Birzan et al. (2008), 

logPcav = (0.35 ± 0.07) logPi4oo + (1.85 ±0.10), (2) 
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where Pcav has units of 10^^ erg s ^ and P1400 has 
units of 10^'* W Hz~^, we estimate the jet power to be 
8.29xl0''3 erg s-\ 

Rafferty et al. (2008) have shown that star forma- 
tion is predominately found in systems where the X- 
ray cooling luminosity is greater than the jet power. 
To test this in A1664, we derived the X-ray luminosity 
by fitting the cooling region with a single temperature 
MEKAL model and using the bolometric flux (0.001- 
100.0 keV) of the model to calculate the luminosity. Sim- 
ilarly, the cooling flow luminosity was calculated from 
the flux of an additional MKCFLOW component with 
the lower temperature set to 0.1 keV. The cooling lumi- 
nosity, 2.45±0.12x lO^"' erg s^^, is calculated as the the 
difference between these two values. The cooling lumi- 
nosity here is three times greater than the estimated jet 
mechanical power foimd from Birzan's relation. 

The uncertainty in the jet power is large. We can check 
our conclusions by using radio power measurements at a 
lower frequency. The radio flux from the nucleus of the 
BCG in the 352 MHz WISH catalogue is 329 mJy (De 
Breuck et al. 2002). With this we can independently 
check the the cavity power estimated using the 1.4 GHz 
flux. Using the 327 MHz relationship of Birzan et al. 
(2008) (where we assume the difference in frequency de- 
pendent flux is negligible), 

logPcav = (0.62 ± 0.08) l0gP327 + (I-H ± 0.17), (3) 

we find a cavity power of 6.67x 10^^ erg s~-^. This value is 
consistent with the cavity power found with the 1.4 GHz 
flux. An additional constraint is found using Birzan's 
relationship between jet power and lobe radio luminosity 
corrected for break frequency to account for aging effects, 

log Pcav = (0.53 ± 0.09) logLrad- 

(0.74 ± 0.26) logi/c + (2.12 ± 0.19). (4) 

Lrad and Pcav are in units of 10^^ erg s~^ and I'c is in 
units of GHz. Correcting for the break frequency reduces 
the scatter in this relationship by ~ 50% compared to 
the monochromatic relations. The VLSS survey gives 
a 74 MHz flux for A1664 of 1.16 Jy. Combining this 
flux measurement with the previous ones, we calculated 
the total radio luminosity by integrating the flux from 
10 MHz to 10 GHz as two broken power law spectra, 
with spectral indices of 0.81 between 10 MHz and 352 
MHz and 1.59 between 352 MHz and 10 GHz. The total 
radio luminosity is found to be 2.01±0.05x 10^^ erg s~^. 
We are using the total flux (lobe plus some core flux), 
so this should be considered an upper limit to the lobe 
luminosity. The shape of the spectrum is defined by just 
three flux measurements. We have therefore taken the 
intermediate flux at 352 MHz to be the rough estimate of 
the break frequency. Based on Birzan et al. (2008), this 
is a reasonable lower limit. Adopting these assumptions, 
we estimate the upper limit on Pcav to be 1.80x10^^ erg 
s~^. Again, this result is consistent with the cooling 
luminosity being greater than the energy output by the 
AGN. 

The scatter in the Birzan power relations is larger 
than the calculated difference between the jet power and 
the cooling luminosity. We cannot exclude the possibil- 
ity that the AGN power exceeds the cooling luminosity. 
However, taking these flgures at face value, the Birzan 



relations give a jet power that is consistently below the 
cooling luminosity. With this, the absence of detected 
cavities, and evidence for cooled gas in gravitational free- 
fall at the center (Wilman et al. 2009), it appears that the 
AGN is currently underpowered compared to the cool- 
ing luminosity. The hot atmosphere is now in a cooling 
phase. 

4.2. Criteria For Cooling Flow Driven Star Formation 

Since A1664 has an unusually blue core in the BCG 
with a star formation rate of 23 M0 yr~^, we have in- 
vestigated whether gas cooling out of the ICM is directly 
feeding the cold gas reservoir and fueling star formation. 
We tested A1664 against the criteria of Rafferty et al. 
(2008) for star formation that require the jet power to be 
less than the X-ray cooling luminosity, a central cooling 
time less than 5x 10^ yr or a minimum entropy threshold 
of 30 keV cm^, and a maximum of 20 kpc separation be- 
tween the optical and X-ray centroids. We have shown in 
previous sections that A1664 has an underpowered jet as 
well as a central cooling time and entropy of 3.5x10^ yr 
and 10.4 keV cm^, respectively. The X-ray and optical 
centroids calculated have a separation of ~ 5 kpc, which 
satisfies all of the criteria. 

There are two other indications we have explored that 
star formation is related to the cooling flow (X-ray cusp) . 
Clusters below the entropy threshold are more likely to 
have high luminosity Ha emission (Cavagnolo et al. 2008; 
Voit et al. 2008), a good tracer of star formation. A1664 
is one the brightest Ha emitters in the ROSAT brightest 
cluster sample (Crawford et al. 1999). Also, the derived 
upper limit on the mass condensation rate (56 Mq yr^^) 
is larger than the star formation rate, although they are 
consistent to within their uncertainties. Rafferty et al. 
(2006) found that the average ratio of mass condensation 
to star formation was approximately four. For A1664, 
the ratio is smaller (^ 2.4). The upper limit on cooling 
being greater than the star formation rate is consistent 
with a cooling flow driven star burst. 

4.3. Gas Deposited by Stripping 

Cold gas and star formation in BCGs are found al- 
most exclusively in clusters with cooling flows (i.e., cen- 
tral cusp of dense gas). Therefore, the origin of this cold 
gas must be linked in some way to the presence of the 
X-ray cusp. Throughout this paper we have argued that 
this link is through condensation from the dense cool 
gas. Another possible link could be through ram pres- 
sure stripping. Because of the high density cusp, the 
cross section to ram pressure stripping is larger in the 
center than in the cluster outskirts. Stripping could then 
be strong enough to drive dense molecular gas out of a 
gas-rich galaxy plunging through the cusp and depositing 
it onto the BCG. 

Blue streams seen in Figure 12 could plausibly be star 
formation in cold gas stripped from plunging galaxies. 
Two of the blue trails appear to end near red ellipticals, 
so their origin is miclear. For the case of the southern 
most trail, Wilman et al. (2006) has suggested that it 
has been dynamically induced by a disturbance from the 
nearby galaxy (point 2 in Figure 2). 

We will now consider the possibility that the 10^° Mq 
cold gas reservoir was built up from gas stripped from spi- 
ral or gas-rich early- type galaxies. Previous observations 
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TABLE 5 

Derived Properties from X-ray & Radio Data 



Sl400 Pl400 S352 P362 Pcav'' Pcav'' Lx LcF r^ool 

(mjy) (1024WHz-l) (mjy) (lO^-* W Hz"!) (10*2 gj.g s"!) (lO^^ erg s"!) (10*2 erg s"!) (lO*^ erg s"!) (kpc) 



36.4±1.2 1.57±0.05 329±9 14.2±0.4 82.9±19.4 66.7±29.8 291±6 46.5±10.8 117 



: 1400 MHz 
: 352 MHz 



have shown that lower density cold gas from a galaxy be- 
gins to be stripped at about 0.6 to 1 Mpc from the center 
of a cluster (Kenney & Koopmann 1999; VoUiikt ct al. 
2001). To determine if any gas will remain in a galaxy 
to be deposited onto the BCG, wc calculated whether a 
face-on, high density disk of molecular gas can survive 
a close passage to the center of the cluster. In order to 
find the radius when stripping should occur, we compare 
the ram pressure stripping force to the maximum gravi- 
tational restoring force of the galaxy parallel to the flow 
as: 



pv 



dz 



(5) 



Here, vj is the radius of the cold gas disk, z is the distance 
perpendicular to the disc, p is the density of the ICM, v 
is the velocity of the galaxy through the ICM, and (Ts is 
the surface density of the cold gas disk. Assuming the 
the galaxy has a NFW potential (Navarro et al. 1997), 
the specific form for our model galaxy is. 
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Fig. 13. — The velocity required at varying radii in A1664 for 
the ram pressure stripping force to ovcrconic the maximum gravi- 
tational restoring force of dense molecular gas in different galaxies. 
From top to bottom, the lines represent typical ULIRG, LIRG, 
large spirals, and small spirals. Along the right side of the plot arc 
the average masses of molecular gas found in the inner dense region 
of the galaxies. Across this range, galaxies towards the bottom are 
easier to strip, but have less gas to contribute to the BCG. Galax- 
ies towards the top are harder to strip and are rare. The dashed 
horizontal line is the velocity dispersion of A1664. 



where r — Vw^ + z^, M is the mass enclosed at r2oo, ?'s 
is the scale radius, and c is the ratio of 7200 to the scale 
radius. For our calculations we assume this ratio to be 
10 to 1. Solving equation 5 for the velocity a galaxy must 
travel to overcome the maximum gravitational restoring 
force per unit area, 
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(7) 



we present in Figure 13 how this varies with cluster ra- 
dius for different galaxy types, assuming the density pro- 
file of A1664 (Figure 9). 

The dotted line in Figure 13 represents a typical small 
spiral galaxy (M = 10" Mq, = 100 Mq pc'^, 
^^200 = 200 kpc, w = 0.3 kpc, z = 1.11 kpc), such as 
IC 342 (Crosthwaite ct al. 2001). Comparing this to the 
cluster's velocity dispersion of 900 km s^^ (Allen et al. 
1995), it is likely the galaxy would lose its gas before 
reaching the center due to the low threshold value, mak- 
ing it nearly impossible for a significant amount of gas to 
be deposited onto the BCG. The solid line represents a 
typical large spiral galaxy (M = 8 x 10^^ Mq, CTs = 400 
Mq pc~^, r2oo = 210 kpc, w = 0.4 kpc, z = 1.37 kpc), 
such as M83 (Crosthwaite et al. 2002) or the Milky Way 
(Scoville fc Sanders 1987). A galaxy of this type would 
be able to hold onto its high density gas until about 80 
kpc from the center. It is plausible then for the gas to 
be deposited onto the BCG. However, according to sim- 
ulations by Benson (2005), in-falling halos in this radial 
velocity range on plunging orbits are rare. On average, 
a spiral galaxy contains only a few 10* Mq of molecular 
gas in the inner 1 kpc. There would have to be hundreds 
of these encounters to build up the amount of cold gas 
seen in A1664. 

Large spiral galaxies can merge and form what is 
known as a luminous infrared galaxy (LIRG). These 
early-type galaxies contain a potential source of cold gas 
that is comparable to the amount of gas found in BCGs. 
The dashed line represents a typical LIRG (M = 3 x 10^^ 
Mq, cTg = 10^ Mq pc~^, r2oo = 250 kpc, = 0.5 kpc, 
z = 1.68 kpc), such as NGC 4194 (Aalto & Huttemeister 
2000). For galaxies in this range, it would be nearly im- 
possible for any of the gas to be stripped, even in the cusp 
of A1664. Similarly, a typical ULIRG (M = 10" Mq, 
fjs = 10^ Mq pc~^, r2oo = 300 kpc, zu=l kpc, z = 2.85 
kpc), represented by the dash-dot line, would also hang 
onto its gas despite ram pressure stripping. Tidal strip- 
ping would be a more viable method of reinoviiig the 
gas from these galaxies. However, tidal stripping has no 
connection to the cooling cusp. If galaxies are depositing 
10^^° Mq of molecular gas in this way, we should see this 
happening in all clusters, not just in cooling flows. 

Based on these calculations, ram pressure stripping is 
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an improbable explanation for the cold gas reservoir in 
A1664's BCG. This leaves us to favor the cooling flow 
driven star formation scenario. 

4.4. Collapsed Cavity System 

About 70 - 75% of cooling flow clusters have detectable 
X-ray cavities (Dunn ct al. 2005). Thus, close to a third 
of cooling flows would either be in a state of pre- or post- 
AGN outburst. When cavities are formed, they sweep 
up low entropy gas into dense shells that may eventually 
break up and sink back to the cluster center. The X- 
ray bar could be the remains of collapsed cavity shells. 
Assuming the mass of the bar to be the minimum mass 
the AGN must had displaced, and further assuming the 
gas density (5.4x10^^ Mq pc^'^) surrounding the bar to 
be the original, undisturbed density of the center, the 
amount of energy needed to displace 3.21x10-^'^ Mq of 
gas would be 3.1x10^^ erg. This energy would corre- 
spond to two cavities 44 kpc in diameter. Energy and 
sizes of this magnitude are common in clusters (McNa- 
mara & Nulsen 2007). 

5. SUMMARY 

We have shown that star formation in Abell 1664 is 

consistent with fueling by the cooling flow. The BCG 
has a blue central color (positive color gradient), a cen- 
tral cooling time below 5x10* yr, a central entropy less 
than 30 keV cm~^, an AGN power less than the X-ray 
cooling luminosity, and a BCG less than 20 kpc from the 



X-ray core. These properties are consistent with other 
star forming cooling flows (Raffcrty et al. 2008). A large 
10^° Mq cold gas reservoir at the cluster center is also 
consistent with a cooling dominated cluster. It is unlikely 
that the gas was deposited by spiral or gas-rich early-type 
galaxies because it would take hundreds of encounters 
for it to be explained by ram pressure stripping alone. 
Other evidence that the cooling flow is directly feeding 
this reservoir comes from the measured mass condensa- 
tion rate of 56 Mq yr~^, which compares well to the 
observed star formation rate of 23 M© yr~^. 

Additional evidence suggesting that A1664 is currently 
between AGN outbursts includes the possible collapsed 
cavity shell remnants we have referred to as the X-ray 
"bar". We have shown that the 3.1x10^^ erg it would 
have taken to displace this gas would have created a pair 
of cavities approximately 44 kpc in diameter, typical of 
what is seen in galaxy clusters. 

We further demonstrated that if gas in the central re- 
gion is cooling at the rates estimated above, we can ac- 
count for the apparent dip in the metallicity found when 
assuming a single temperature for the medium. 
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